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Large-scale inhomogeneities may improve the cosmic concordance of supernovae 
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We reanalyze the supernovae data from the Union Compilation including the weak lensing effects 
caused by inhomogeneities. We compute the lensing probability distribution function for each back- 
ground solution described by the parameters Qm, Qa and w in the presence of inhomogeneities, 
approximately modeled with a single-mass population of halos. We then perform a likelihood anal- 
ysis in the space of FLRW-parameters and compare our results with the standard approach. We 
find that the inclusion of lensing can move the best-fit model significantly towards the cosmic con- 
cordance of the flat ACDM model, improving the agreement with the constraints coming from the 
cosmic microwave background and baryon acoustic oscillations. 
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Introduction. In the standard approach super- 
novae (SNe) observations are analyzed in the frame- 
work of homogenous Friedmann-Lemaitre-Robertson- 
Walker (FLRW) models. However, the universe is known 
to be inhomogenous, showing a distribution of large 
galaxy clusters and filamentary structures surrounding 
much emptier voids of size sw 10 — 100 Mpc. A known 
effect of these structures on any set of standard candles 
is weak gravitational lensing [1 . Weak lensing can cause 
either brightening or dimming of the source depending 
on whether the matter column density along the line of 
sight is larger or smaller than the FLRW value. 

The fundamental quantity describing this statistical 
magnification is the lensing probability distribution func- 
tion (PDF). The lensing PDF is specific both to the given 
FLRW model, and to the particular spectrum of inho- 
mogeneities introduced. It is not currently possible to 
extract the lensing PDF from the observational data and 
we have to resort to theoretical models. Two possible 
alternatives have been followed in the literature. A first 
approach (e.g. Ref. [5]) relates a "universal" form of the 
lensing PDF to the variance of the convergence, which 
in turn is fixed by the amplitude of the power spectrum, 
a§- Moreover the coefficients of the proposed PDF are 
trained on some specific N-body simulations. A second 
approach (e.g. Ref. [5]) is to build a model for the inhomo- 
geneous universe and directly compute the relative lens- 
ing PDF, usually through time-consuming ray-tracing 
techniques. The flexibility of this method is therefore 
penalized by the increased computational time. 

In this work we use another approach, based on the 
stochastic modelling of the inhomogeneities introduced 
in Ref. [¥| . This method combines the flexibility in mod- 
elling with a fast performance in obtaining the lensing 
PDF. To compute one lensing PDF, the numerical im- 
plementation turboGL 0.4 (4] takes, with an ordinary 
desktop computer, a time of order of a second. This 



speed performance makes it feasible to do an ab initio 
likelihood analysis in the space of FLRW- models endowed 
with inhomogeneities. In this letter we will perform such 
an analysis for the Union SNe Compilation [5]. 

Setup. We will treat inhomogeneities as perturba- 
tions over the FLRW model which is parametrized as 
usual by the present Hubble expansion rate H = lOO/i 
km s _1 Mpc -1 , the present matter density parameter Q,m 
and the present dark energy density parameter fl\ and a 
constant equation of state w. We fix the radiation den- 
sity to fifl = 4.2 ■ 10~ 5 h~ 2 . For inhomogeneities we use a 
"meatball" model [6 consisting of randomly placed spher- 
ical halos made of ordinary and dark matter. In princi- 
ple these halos need not be virialized, and the spherical 
symmetry assumption is not very restrictive. As was ex- 
plained in Ref. , the weak lensing properties of a given 
universe model can be described by a set of matter dis- 
tribution projections (z-dependent column densities) on 
a small number of independent redshift slices. For such 
projections any local density contrast, such as a long fil- 
ament seen edge on, looks roughly like a spherical halo. 

Here we use a simple single-mass halo model which 
is completely parametrized by the comoving distance be- 
tween halos A c , the halo proper radius R p and the density 
profile. We choose the latter to be the Navarro-Frenk- 
White profile [7] with a concentration parameter c ~ 6.7 
and we assume that the halos have virialized with a con- 
trast of 200 at a redshift z v ; r , whereby (for a given z v - n ) 
the corresponding R p can be taken constant. The halo 
mass is related to the comoving density n c = A~ 3 by 
Pc^Im = M n c . For numerical values we explored the 
range A c = (5.4, 9.0, 12.6) /i _1 Mpc and correspondingly 
M = (0.44, 2.0, 5.6) IO^-^mMq for z vir = 0.8, and 
z vir = (0, 0.8, 1.6) for A c = 12.6 h~ l Mpc. The numerical 
value of R p depends on the background matter density 
at z V j r . For the ACDM model the previous range of z v i r 
corresponds to R p ~ (0.9, 0.7, 0.5) ft _1 Mpc. 
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Lensing. The meatball model incorporates quantita- 
tively the crucial feature that photons can travel through 
voids and miss the localized overdensities. This feature is 
not present, for example, in "swiss-cheese" models where 
the bubble boundaries are designed to have compensat- 
ing overdensities. Such models have indeed been shown 
to have on average little lensing effects [HIE]- The key 
quantity in all our analysis is the lens convergence k, 
which in the weak-lensing approximation is given by 
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a(t) and t(r) correspond to the FLRW model, r s — r(z s ) 
is the comoving position of the source at redshift z s 
and the integral is evaluated along the unperturbed light 
path. Also, fk{r) — sin(rv / fc)/v / fc, r, sinh(r-\/— k) j\f—h 
depending on the curvature k >, =, < 0, respectively. 

Neglecting the second-order contribution of the shear, 
the shift in the distance modulus caused by lensing is 
expressed solely in terms of k: 



Am(z) = 51og 10 (l - k(z)) . 



(2) 



Equations (T][2) show that for a lower-than-FLRW col- 
umn density the light is demagnified (e.g., empty beam 
5 = —1), while in the opposite case it is magnified. 

In Ref. jl] a fast and easy way to obtain the conver- 
gence PDF for these meatball models was derived. In 
short, the formula for the convergence Eq. ([T| is replaced 
by a discretized probabilistic expression: 
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Here Ku m is the convergence due to one halo, at 
a comoving distance r^, which the photon path in- 
tercepts with an impact parameter b m , n\i m = 
G{r u r s ) j^%m where Pi (x) is the local 

halo density and pM is the FLRW matter density. In 
practice one divides the comoving distance r s to the 
source and the radius R of the halo into bins of widths 
R <C Ari r s and Ab m <C R and lets the centers of 
these bins define the allowed values for r and b. The 
quantity ki m in Eq. pj is a Poisson random variable 
of parameter ANi rn = n c AVi rni which gives the ex- 
pected number of halos within the phase space volume 
AVi m — 27r6 m A6 m Ar i . That is, Eq. (|3| defines a conver- 
gence as a function of a configuration {ki m } of halos along 
an arbitrary line of sight from the observer to the source. 
The lensing PDF in the distance modulus P w \(Am, z s ) 
is then constructed from a large sample of random con- 
figurations {k im } using Eqs. ^ and Note that the 
expected convergence computed from Eq. ^ is zero, con- 
sistent with photon conservation in weak lensing, because 




Figure 1: PDFs for a SN with a = 0.25 mag at z 3 = 1.5 in 
the ACDM model endowed with halos specified by z v i r = 0.8 
and A c = 12.6 h~ 1 Mpc. The dotted histogram represents the 
lensing PDF, the dashed line the SN PDF and the solid line 
the full likelihood. The dotted curve is described in the text. 



for a Poisson distributed variable the expected value co- 
incides with its parameter. 

Likelihood function. After the raw lensing PDF 
P w i(Am) has been computed for a given set of FLRW- 
parameters and redshifts, it still has to be convolved 
with the intrinsic source brightness distribution Pi n : 
P(Am,z s ) = f dyP w i(y,z s )P in (Am-y). We take P in to 
be a gaussian in the distance moduli. The actual intrin- 
sic distribution should be a universal function if the SN 
are similar at all distances. However, following Ref. [5], 
we will combine all observational (gaussian by assump- 
tion) uncertainties in quadrature with the intrinsic dis- 
tribution, whereby P lVL becomes an effective distribution 
specific for each SN event Pi n (^) ~> Psn(x, (Jj). The 
likelihood function for a single SN-observation is then 

Li{^) = JdyPw\(y,z i )P S N(Am i -^-y,a i ), (4) 

where Ami = "Vi — ^o,i is the observed magnitude 
and the corresponding FLRW prediction is related to the 
luminosity distance by m t i — 51og 10 (^(z^/lOpc. 
The parameter fi is an unknown offset sum of the SNe 
absolute magnitudes, of fc-corrections and other possi- 
ble systematics. Note also that Li inherits the vanishing 
mean of P w \ and that its variance is simply given by the 
sum of the variances of the convolving PDFs. 

We define the total likelihood function as the product 
of all independent likelihood functions in the data sam- 
ple, further marginalized over \i: 



L(n M ,^A,w) 



d/iIliLi(iJ,) . 



(3) 



Since /x is degenerate with log 10 Hq we are effectively 
marginalizing also over the expansion rate of the uni- 
verse. A replacement of P w \(y,z) by a cosmology- 
independent gaussian with a variance [5] a = 0.093z, 
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Figure 2: 1,2 and 3cr confidence level contours on w and Qai", 
for a flat wCDM universe with halos specified by z v [ r — 0.8 
and A c = 12.6ft~ 1 Mpc. The results using the full likelihood 
of Eq. (|5f are shown as filled contours and the best-fit model 
with a circle. The results using the gaussian Gi are shown as 
dotted lines with a triangle, the ones using the unlensed Psn 
are shown as dashed lines with a square and correspond to 
the ones of Ref. [5] (without systematics). 



would reduce Eq. (|5J) to the form used in the analysis of 
Ref. [5]. Typical forms of P w j, Psn an d Li(/j, — 0) have 
been illustrated in Fig. [I] Also shown for later use is Gi, 
which is a gaussian with the same variance of Li(0). 

Results. We run a global likelihood analysis using the 
formula (|5| for two different setups: first in the (Qm , w)- 
space for flat (fifc = 0) wCDM models and second in the 
(Om, f^A)-space for a non-flat ACDM model (w = — 1) 
using the Union SNe Compilation of Ref. [5 . We show 
our results in Figs. [2] and [3] as confidence level con- 
tours for x 2 = ~~ 2 log L. For comparison we have per- 
formed the analysis also using the standard Psn dis- 
tribution (as done in Ref. [5]) and the distribution Gi. 
The idea for using Gi is that it takes into account the 
cosmology-dependent extra dispersion coming from lens- 
ing, but neglects the skewness of the true distribution. 
So, the contours relative to Gi give an idea of how much 
of the difference from the standard analysis comes from 
the widening of the intrinsic distribution, and how much 
from the skewness of the actual PDF. As it is evident 
from Figs.[2}|3] the la contours are basically determined 
by the cosmology-dependent widening, whereas the skew- 
ness starts to be relevant between the 2 and 3a levels. We 
point out that our results are essentially unaffected if we 
add a constat a sys to the Oi of Eq. Q in order to have 
the same reduced x 2 using Gi and Psn- 

Our most important result, clearly evident from 



Figure 3: 1, 2 and 3<r confidence level contours on Qm and J1a 
(i.e., allowing for non-zero curvature) for ACDM (w = — 1) 
with halos as in Fig. [2] Note that, as also in the previous plot, 
the new best-fit points lie on the la confidence level contour 
relative to Psn- Labeling as in Fig. [2] 



Figs. [2f|3j is that the inclusion of lensing effects in 
the likelihood analysis significantly moves the best-fit 
model, from {n\ IlW v ) = (0.38,-1.4) and {9^,0%) = 
(0.41,0.94), towards the cosmic concordance of the flat 
ACDM model, therefore improving the agreement with 
the constraints coming from cosmic microwave back- 
ground (CMB) and baryon acoustic oscillations (BAO). 

To further explore this behavior we studied how the 
new best-fit model positions (Q* M ,w*) and (Q* M ,fl A ) de- 
pend on the halos mass M and the virialization redshift 
z v j r . For a fixed halo mass, higher values of z v ; r give 
denser halos with smaller radius R p and higher lensing 
corrections to the likelihood. As explained before, the 
numerical value of R p depends on the background model 
and we use the ACDM as a reference model to convert 
z v i T into R p . A fit for A c = 12.6/i _1 Mpc then gives: 

{n* u ,w*) = - 1.5e- 2 - 3R "],w v [l - 0.94e- 1 ' 7 ^]) 

(^,^) = (n^[l-0.8e- 1 ' 2 ^] ) nX[l-0.54e- 1 ' 1 ^]) 

(6) 

where R p > 0.5 is in units of h^ 1 Mpc. If we fix z v - lr , 
higher values of M (or equivalently A c ) give a universe 
made of larger clumps with larger voids giving therefore 
stronger lensing corrections. A fit for z vil = 0.8 gives: 

(n* M , w*) = - 0.13 M - 47 , w v + 0.45M - 33 ) 
(Sl* M , Q* A ) = (fi& - 0.15 M - 29 , n\ - 0.25 M - 26 ) (7) 

where M < 1 is in units of 5.6 • 10 u h- 1 n M M . 
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The general trend favoring models with smaller £1 m fol- 
lows from the fact that lensing effects in general make the 
fit slightly worse with than without lensing [10]. The ef- 
fect comes both from the cosmology-dependent widening 
and from the skewness of the distributions, and it is obvi- 
ously more pronounced for larger matter densities. This 
can be seen directly from Eq. ([TJ, where the magnitude 
of the lensing effects is explicitly seen to be proportional 
to f2j\/. The overall movement of the best-fit model then 
follows the degeneracy of the FLRW models. 

Discussion. Our halo model was designed to capture 
the most important effects of the weak gravitational lens- 
ing by the nonlinear large-scale structures. In particular 
the large voids that dominate the late-time universe were 
imposed in the model by concentrating all matter into ha- 
los. Accordingly, we chose the halos to have the mass of 
a very large cluster, i.e. of order 10 li h~ 1 M & , which then 
corresponds to an interhalo distance of order 10/i _1 Mpc. 

Given this result, it is natural to ask if our toy model 
could also give a reasonable approximation to the ob- 
served power spectrum. This is not entirely obvious, 
because weak lensing and power spectrum probe some- 
what different aspects of the inhomogeneities. We found 
that our single-mass halo model tends to concentrate too 
much power onto the interhalo distance scale, when com- 
pared to the non-linear correction to the ACDM spec- 
trum provided for example by the halo model of Ref. [IT] . 

It will clearly be interesting to improve the modelling of 
the power spectrum by adopting a more realistic halo dis- 
tribution function /(M, z), and we plan to pursue this in 
future work. However, this paper was devoted to explore 
the extent to which lensing can change the supernovae 
results and this is best done by adopting the simplest 
single-halo model with its few parameters. In any case, 
the choice of the mass function /(M, z) is not that simple; 
even after fitting the power spectrum well, an efficient 
lensing requires modelling the voids and filaments that 
are described by higher order correlation terms. This can 
in principle be done in the current approach by introduc- 
ing additional large scale variations to the background 
density from which standard halo functions are drawn. 

Finally, given a large enough SNe dataset one could 
in principle measure the lensing PDF. However, to do 
this properly one would have to understand the selection 
effects that could for example cut the high magnification 
tail of the PDF, sizably biasing the average convergence, 
variance and skewness. 

Conclusions. We have presented a reanalysis of the 
supernovae data from the Union Compilation including 
the lensing effects caused by inhomogeneities. Unlike 
in the analysis of Ref. |5j [12] , where the lensing effects 
are accounted for by adding in quadrature a small z- 
dependent variance to the other statistical and system- 
atic errors, we compute the actual probability distri- 



bution functions for each different FLRW-model with a 
spectrum of inhomogeneities designed to mimic the ob- 
served large-scale structures. In particular, large voids 
that dominate the late-time universe are imposed on the 
model by concentrating all matter into halos of mass 
of order 10 14 h~ 1 M & . We found that including inho- 
mogeneities significantly changes the likelihood contours 
(the likelihood peaks, for instance, move of around la) 
and clearly improves the concordance of the supernova 
data with the CMB and the BAO, which may be used to 
strengthen the case for the standard ACDM model. 

One should be reminded that our findings could change 
if other effects caused by large-scale inhomogeneities are 
introduced, e.g., selection or redshift effects. It also re- 
mains to be seen how a more realistic inhomogeneous 
distribution, providing a better fit to the matter power 
spectrum, would affect these weak- lensing corrections to 
the SNe contours. 
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